UR-1298 
ER-40685-747 
hep-ph 930120S 
March 1993 



Spherical Quantum Chromodynamics of Heavy Quark Systems 



K. S. Gupta, S. Guruswamy and S. G. Rajeev 



Department of Physics and Astronomy 
University of Rochester 
Rochester, N.Y. 14627 
e-mail: gupta@uorhep, guruswamy @uorhep and rajeev@uorhep 



Abstract 

We propose a model for Quantum Chromodynamics, obtained by ignoring the angular 
dependence of the gluon fields, which could qualitatively describe systems containing one 
heavy quark. This leads to a two dimensional gauge theory which has chiral symmetry 
and heavy quark symmetry. We show that in a light cone formalism, the Hamiltonian of 
this spherical QCD can be expressed entirely in terms of color singlet variables. Further- 
more, in the large N c limit, it tends to a classical hadron theory. We derive an integral 
equation for the masses and wavefunctions of a heavy meson. This can be interpreted as a 
relativistic potential model. The integral equation is scale invariant, but renormalization 
of the coupling constant generates a scale. We compute the approximate beta function of 
the coupling constant, which has an ultraviolet stable fixed point at the origin. 
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1. Introduction 

Two dimensional Quantum Chromodynamics (QCD) is fairly well understood in the 
■jtj- expansion, following the work of 't Hooft, Witten and others [1], [2], [3]. Furthermore, 
it can be solved by numerically diagonalizing the hamiltonian [4] . The two methods are in 
satisfactory agreement with each other. In this paper we will study an approximation to 
four dimensional QCD in which the gluon field is independent of two angular variables, so 
that two dimensional methods can be applied. Although this could be studied just as a toy 
model for QCD, it might also be a model for hadrons containing one heavy quark. These 
systems have attracted much attention recently with the discovery of a new heavy quark 
symmetry [5], [6], [7]. Such a hadron is similar to an atom, the heavy quark being like 
the nucleus with the light quarks orbiting like electrons around it. The main difference is 
that the light quarks have to be described relativistically. Still, we should be able to treat 
these systems in the spirit of the Hartree-Fock or Thomas-Fermi methods [8]. The gluon 
field produced by the heavy quark is spherically symmetric. That of the light quarks is not 
spherically symmetric because the light quark wavefunctions depend on angles. Still, we 
should be able to do a 'spherical averaging' as in Hartree-Fock theory and approximate the 
current density of the light quarks by their average over the angular variables. The average 
of the angular components of the current density will vanish. This spherical averaging is 
equivalent to putting Aq = = and assuming that the remaining components A are 
independent of the angles. This will lead to an effectively two dimensional theory ( although 
without translation invariance) in light cone coordinates. This 'Spherical QCD' respects 
heavy quark symmetry, chiral symmetry of light quarks and scale invariance. (The scale 
invariance should be broken by quantum effects.) It should be possible to solve this theory 
numerically, using methods similar to those used for two dimensional QCD [4]. This is 
vastly simpler than solving the full 3 + 1 dimensional problem numerically. In this paper, 
we will study this spherical QCD analytically using the expansion. 

In this paper, we will view Spherical QCD as a model to study theoretical aspects of 
QCD such as the expansion and asymptotic freedom. We have a model that is solvable 
in the large N c limit, in which the beta function can be calculated in a non-perturbative 
approximation. Whether our approximate theory can make reliable numerical predictions 
can only be seen from future work. 

't Hooft's original approach to the expansion involved summing planar Feynman 
diagrams. This method is difficult to generalize to our situation. More recently, another 
method has been used to construct a bilocal quantum field theory of hadrons equivalent to 
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2DQCD for all N c [9] . The parameter plays the role of h in this hadron theory ( N c must 
be an integer for topological reasons) so that the large iV c limit of QCD corresponds to a 
classical theory of hadrons. The field equations of this theory are certain non-linear integral 
equations: the large N c limit of QCD is not a free field theory. Semiclassical expansion 
around the classical solutions will describe the hadron masses and wavefunctions in the 
expansion. The most obvious static solution is the vacuum, and small fluctuations around 
it (mesons) are described by a linear integral equation, in agreement with 't Hooft. This 
point of view also allows for baryons, which are then static solutions that deviate by a 
large amount from the vacuum [9]. 

This approach to two dimensional QCD can be generalized to Spherical QCD. One 
of our results is a linear integral equation for the heavy-light meson masses and wave 
functions. We can understand our result ( just like 't Hooft's integral equation) as a 
relativistic potential model. Ordinarily, such a one particle point of view would not be 
allowed in this highly relativistic situation due to effects of virtual light qq pairs. In the 
meson picture, these are the effects of virtual light mesons. However, in the large iV c limit 
such effects are suppressed, so that a one particle description is allowed. Thus we will get 
a light quark moving in the field of the heavy quark. The Dirac equation can be separated 
in light cone coordinates, reducing the problem to a linear integral equation. (That light 
cone coordinates are better suited to nonperturbative problems in QCD have also been 
suggested for other reasons [10], [11]. ) It should be possible to fit meson masses to a 
potential with a small number of parameters as a phenomenological test of our approach. 

In the next section (Sec. 2 ) we will derive the action of spherical QCD; it will be 
convenient to use a light cone coordinate system. In Sec. 3 we will show that this theory 
can be written entirely in terms of a set of color singlet fields representing hadrons. The 
commutation relations and equations of motion of these variables are derived. In Sec. 4 
we show that in the large N c limit, this tends to a classical theory whose Poisson brackets, 
constraints and equations of motion are derived. By expanding around the vacuum solution 
we get a linear integral equation for the masses and wavefunction of heavy mesons. In Sec 
5. we have analysed the ground state of the wave equation in a variational approach. A 
renormalization of the coupling constant is necessary to keep the ground state energy finite 
as the cut-off is removed. We compute the approximate beta function and show that it 
has an ultraviolet stable fixed point at the origin. Finally, in the appendix we describe 
how the constraints on the color singlet variables can be derived, in a simplified context. 
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2. QCD in light cone coordinates 

The action of four dimensional QCD with N c colors, n / massless quarks and Nf heavy 
quarks of mass M is 



where and are N c x N c hermitian matrices. Here, the color indices are suppressed 
and a = 1, • • • nf, A = 1, • • • Nf. We are interested in systems that contain only one heavy 
quark, which can be assumed to be at rest at the origin. Conventional spherical polar 
coordinates t, r, 9, (p are cumbersome in this context due to the negative energy states of 
the light fermions and the dependence of the Dirac sea on the gluon fields. We find it 
most convenient to introduce a sort of light cone coordinate system u,r,6,(j) centered at 
the heavy quark: 



Note that this is not an orthogonal coordinate system; yet the Dirac equation is separable 
in this system. The surface u =constant is a past light cone with apex at the origin. We 
will use a canonical formalism in which initial data is given on this null surface. The vector 
field is in fact just the usual time translation, so that its conjugate variable is energy. 
But for fixed u is not the same as the radial vector in spherical polar coordinates. 
Instead, it is a null vector pointing outward along the surface u =constant: 



We will expand all the fields in appropriate spherical harmonics in the angular vari- 
ables. The essential approximation we will make is to include only the s-wave states of 
the gluon fields. This amounts to replacing the current density of the light quarks by its 
average over the angular variables, so that the mean gluon field they produce is spherically 
symmetric. ( This concept of a 'spherical average' was introduced into atomic physics by 
Hartree-Fock and is extensively used in calculations of atomic energy levels and wavefunc- 
tions [8].) We will be able to eliminate the gluon fields and half the degrees of freedom of 



S 




u = t + r, ds 2 = du(du - 2dr) - r 2 (d6 2 + sin 2 9d<p 2 ). 



(1) 
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the light quarks, to get a theory of quarks interacting through a selfconsistent potential. 
The resulting nonlocal lagrangian will have chiral symmetry, heavy quark symmetry and 
scale invariance, and defines our 'Spherical QCD'. 



Let us first simplify the gluon Lagrangian. After spherical averaging, the angular 
components of a spin one field must vanish, and A u and A r would be independent of the 
angles. Furthermore we can choose the gauge condition 

A r = 0, (3) 

leaving only the component A u = A. In this case the Yang-Mills field strength is 

F^dx 11 A dx v = d r A u dr A du (4) 

and the Yang-Mills action is 

S YM = - I tr[0 r A] 2 r 2 drdu. (5) 
a J 

2 

( Here a = |-.) In this gauge there are no propagating components for the Yang-Mills 
fields since there is no derivative of A with respect to u in the action. 

Let us now turn to the light quark action. The free Dirac operator in the usual 
spherical polar coordinates is [12] 

_ . , d , d 1 K ^ 
l7 .0 = ?7o _ + a {_ + } 6 

ot or r r 

where 

K = S-L + 1, a 2 r = l [a r ,K}+ = 0. (7) 

As usual S are the Dirac spin matrices, a r = 7o7r, and L is the orbital angular momentum 
operator. In our coordinates, 

• r/- ^ d , 8 1 K 

i 1 -d = i l0 [{l + a r )— + a r {— + }. 8 

au or r r 

It will be convenient to diagonalize a r and 7 5 (we define 7 5 such that 7 | = 1): 

75^fc,e,M = €V k,e,LL a r Vk,e,n = I^Vk,e,n Kv k ,e,fi = kv k ,e,-fi- (9) 

The eigenvalue k of K is a positive integer. The spinorial harmonics Vk, 6 ,n are orthogonal 
and are normalized such that 

J dn\v k , e J 2 =4vr. (10) 



Then ip can be expanded in this basis (suppressing color and flavor indices on ip and x): 

V> = Yl /( 8n ) r Xk,eA r ) V k,e,»- (11) 

/n=±l,fc>0,e=±l ^ ^ ' 

The factor ^^tt) r * s c h° sen so as to simplify later expressions. In this basis, a r = 03,-K" = 
ko\ and —ia r K = /ccr 2 - The free Dirac action is then 

/l f t d d k 

#7 • d]^d 4 x = - 2^ / drduxi e [*(l + + + -^Xfce, (12) 

where » . = f ^ fc,e ' 1 ) . Including the gauge fieid amounts to replacing r X by r X - iA: 

\Xk,e,-lJ 

$[17 ■ d + J ■ AtycPx = -J2 J drduxlM 1 + ^)[-^ ~ iA] + ta 3 — + -a 2 ] X ke. (13) 

fc,e 

The heavy quark action can be simplified by the transformation [5], [6] 

* -> e iMv ' K #, (14) 

so that in the limit M — > oo 

y # A [ry • d + 7 • A + M]^ A d 4 x -> i ^ • [0 - L4]^ A ci 4 x. (15) 
Here v is the heavy quark velocity. Now we put 

** = 7(^- QAi (16) 

the spinor Q satisfies 7 ■ vQ = Q [5], [6] and depends only on r. In our coordinate system, 
v = (1, 0, 0, 0) so that this becomes just 

ijQ^-iAWrtu. (17) 

(We may suppress flavor indices as well when they are not essential) . 

Thus the complete action is 

S =i J Q^[d u - iA]Qdrdu tr[d r A] 2 r 2 drdu 



In this action, A and Xk,e,-i do not have derivatives with respect to u. Hence they 
do not propagate and can be eliminated by solving their equations of motion. We have, 

dr^drAa] = Ct Pci (18) 

where, 

Pa = Y^ Xfc, e ,l*aXfc,e,l + Q^t a Q, (19) 

k,e 

and t a are the color matrices. We normalize them such that 

trtatp = ^ a p. (20) 

Half the light quark degrees of freedom also do not propagate: 

d k 

-ifrXk,e,-l +i~Xk,e,i =0. (21) 



Then, 



where 



and 



/■OO 

A(r) = a G(r,r')p(r')dr' (22) 
Jo 



G(r,r') = -min(i 1) (23) 



Xk,e,-i = kd- 1 ^ = f ^ Xk ,eAr')dr'. (24) 

r J oo r 

The boundary condition are that as r — > oo the fields vanish. That is, there is no incoming 
radiation from past null infinity. 

We can now eliminate A, Xke,-i to get an action depending only on Xkei and Q: 
S =i Qdrdu + / drduxl el Xkei 

ke 

+ o XI / dr duxl £l [id r + k 2 r~ 1 (id r )~ 1 r~ 1 ]xkei ~ j dr'drdup a (r)p a (r')G(r,r') 
ke J J 

: J dudrdr'G(r,r')Q^(r)t a p a (r')Q(r) + a J dudrdr' G(r 7 r')x\ el (r)t a p a (r')xkei(r) 



ke 

+ a 



This action defines Spherical QCD. It is manifestly invariant under U{nf)i, x U(nf) R ( 
chiral symmetry), SU(2N — f) ( spin and flavor symmetry of heavy quarks) and scale 



transformations ( the only coupling constant g 2 is dimensionless) . The first two terms 
simply say that Q,Xkei are canonically conjugate to Q,x\ e i- 

[Qt(r), Q(r')]+ = S(r - r'), [ X \ el (r), Xfcvi(r')] + = 5 kk ,5 £e ,5(r - r'). (25) 



The remaining terms determine the hamiltonian: 
H = ^2 J dr xt el hXkei + J drdr' ' p a (r)p a (r')G(r,r') 

ke 

a J drdr'G(r, r')($ (r)t a p a (r')Q(r) - a J drdr'G(r, r')xt el (r)t a p a (r')xkei{r), 



where 



h = --(id r + k^-^idr)-^- 1 ) 



is the single particle hamiltonian of the light quarks. Since the hamiltonian H is conjugate 
to u, it is just the energy. The above anticommutation relations along with the hamiltonian 
define Spherical QCD. 

Apart from the fact that h and G have different explicit forms ( and the presence of 
the heavy quarks), this is exactly like the hamiltonian of 2DQCD, after eliminating the 
gluons and half the quark fields. Thus from this point, we can follow the method described 
in reference [9] to turn this into a theory of bilocal color singlet ( meson) fields. 

3. Hadron Theory 

Let us define the color singlet bilinears, 

H i(r',r) = ±-tf I (r') X i(r) , fljV, r) = J^VjQjM, 

Pj{r,r') = ±-Q^(r)Qj(r>) , M;(r,r>) = ±- : X A {r)Xi{r') ' 

where % = (k, e, a) ranges over spin,chirality and flavor of the light quark while I labels 
the spin and flavor of the heavy quark. P is normal ordered with respect to the trivial 
vacuum; after the transformations we made the heavy quark has no negative energy states. 
The normal ordering of M is defined with respect to the vacuum obtained by filling the 
negative energy states of the one-particle hamiltonian h = — \[id r + A; 2 r~ 1 (z9 r )~ 1 r~ 1 ]. In 
fact 

± : X l hr) Xj (r') := i/ (r ) X](r ') + l _ [e{r y) _ 5 ( r , r')]<*j 



*c 1 'c 

1 r A 1 
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Here e(r, r') is the operator that is +1 on the eigenstates of h of positive energy and —1 
on those of negative energy: 



poo M 

s(r,r') = / d\u x (r)ul(r') - / d\u x (r)u* x (r'). 

JO J -oo 



(26) 



The eigenstates u\oi h are defined as below: 



1 



-[id r + k r (id r ) r ]u\(r) = Xu\(r) 



/oo 
u x {r)u* x {r')d\ =8{r-r'). 
-oo 

Also, we will often use the projection operators to the positive and negative energy states 

1 r°° 

S+(r,r') = -[S(r-r') + e(r,r')}= u x (r)u* x (r')d\ (27) 

2 Jo 



5-(r,r') = -[S(r — r') — e(r, r')] = / u x (r)u* x (r')d\ 



(28) 



We will rewrite the commutation relations and the hamiltonian entirely in terms of 
these color singlet variables. The commutation relations are 



[H{( r y),H${ a ,j 

[Hl(rS),M>(8,J 

[H?(r,r'),P£(s,s' 
[H$(ry),Ml(s,s' 

[p/(ry),pf( a y 

[Mi(r,r'),Mt(s,s' 



1 



^-[Pj(r, s')5l5{r' - s) - M({s, r^s' - r) 



1 



ff/( S ,r')<^(r- S '), 



1 



= irl p £(r, s')5?5{r' - s) - If (s, r')5 J L 5(s> - r)], 



K( „ J\xJ X ( J 



1 

iv; 



[M/(r, s')5?5(s - r') - M« (s, r')^5(r - s'), 



+ -(-e(r, s')<J(s - r') + e(s, r')S(r - s'))SfSf}. 



The c-number terms proportional to e(r, r') and 5_ arises because of the normal ordering 
of M and P. 
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The hamiltonian can be expressed in terms of these bilinears after a Fierz reordering 
using the identity for the color matrices: 



L aq l as 2 q ' 



(29) 



(We are using the gauge group U(N C ) rather than SU(N C ), as is usual in large iV c ap- 
proaches [1].) We get, with a = aN c 



It should be possible to find the eigenvalues of this hamiltonian by numerically diag- 
onalizing a finite dimensional approximation to it. This is the method used in Ref. [4] to 
solve two dimensional QCD. 

Now we note that N c appears only as an overall constant in the commutation relations 
and the hamiltonian. In particular, the equations of motion are independent of N c . This 
means that the limit N c — > oo is a classical limit in which the commutators get replaced 
by Poisson brackets. 

The large N c limit of spherical QCD is a classical hadron theory. The semiclassical 
expansion around a solution of this classical theory will be equivalent to the expansion 
of spherical QCD. 

4- Classical Hadron Theory 

The Poisson brackets obtained by the correspondence principle 



H 




[A,B] = ±{A,B} 
10 



(30) 



where plays the role of %. Thus, 

i{H^r,r'),H^(s,s')} = Pj(r, s')5{5{r' — s) — M?(s, r')6j5(s' — r) 

- 6-(8y)S 3 i S I jS(8' -r), 
i{Hl (r, r>),Ml(s, s')} = H{(r, s')^' - s), 
i{H!(r, r'),P£(s, s')} = -H?(s, r')<^(r - s'), 
i{H${ry),P J K {s,s')} = H$(r,s')Sj6(r' - s), 
i{H?(r,r>),Ml(s,s>)} = —H 3 ^(s, r')8\.8{r — s'), 
i{P/(r, r'), Pf ( S , *')} = P/(r, S ')5f 5(r' - a ) - Pf ( s , r')<^(s' - r), 
i{M/(r, r'), M, fc (s, *')} = Mf (r, a 7 )*?^ - r') - M z k (s, r')5j5(r - s'), 

+ ±(-e(r, s')<^ - r') + e(s, r')5(r - J))S}6}, 

and the hamiltonian 

H = J h(r, r')Mi(r, r')drdr' + ^a J G(r, r')[Pj(r, r')Pj \r' » 

+ #/ (r, r')^(r', r) + if^(r, r')#/ (/, r) + M](r, r')M/(r', r)]drdr' 
define a classical theory. The equations motion of this classical theory can be worked out 
by a straight forward ( but tedious) calculation of the Poisson brackets of the observables 
with the hamiltonian. We get 

i ^d~^ = l dc J ds ' H t( r ' s')G(s',r)5-(s',r') + J ds'H^r, s')h(s',r') 
-\aj ds'G(r',s')[Mi(s',r')Hi(r,s') + pi(r,s')H?(s\r'K 

+ ±aj ds'G{s\r)[Pjy S ')H J i {s'y) + Mi{s'y)H I j {r,s% 

dMUr r') 1 r 
i = -a ds'[G(s', r')Mi(s', r')e(r, s') - G(r, s')Afj(r, s')e(s' , r')] 

-^aj ds'G(s'y)[M l j (s?y)Mi(r, s') + Hj(s', r')H\\r, s')} 
+ \aj ds'Gys')[M]{s\r')Miys')+H^s\r')HYys')] 
+ J ds[Mi(r, s)h( S y)-h(r, s)M}{s,r% 
^"du^ = ¥ I ds'G(r,s')[Pk(r,s')P^(s'y) + Hl(r 1 s')Hl\s' 1 r')] 



\a J ds'G{s\r')[Pk{r,s')P«{s\r') + Hl{r,s')H^s'y)\. 
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These equations of motion however do not describe the classical hadron theory com- 
pletely. The bilinear s satisfy some constraints, whose origin is described in the appendix 
in a simplified fermion theory. The constraint for M has been derived in the appendix. It 
becomes in our present notation, 

J Mj(r, r')M J k (r', r")dr' - 1 J [e(r, r')M l k {r\ r") + M l k (r, r')e{r' , r")]dr' 

+ J H} l (r,r')H[(r',r")dr' = 0. 

Similar arguments show that 

J H*(r, r')M}(r', r")dr'- J Hj (r, r')6+(r', r")dr 

+ J Pj(r, r')Hl{r', r")dr' = 

and 

J Hl{r,r')H]\r',r")dr' + J P^(r,r')P? \r' \r")dr' - Pj(r,r") = (31) 

The equations of motion and the constraints above define the classical hadron theory 
equivalent to spherical QCD in the large iV c limit. 

Any static solution to these equations of motion and constraints can be used as the 
starting point for a semiclassical expansion of spherical QCD. The most obvious solution 
is the vacuum: 

Mj(r, r') = #/(r, r') = Pj(r, r') = 0. (32) 

The small fluctuations around this vacuum will describe the mesons in our approxi- 
mation. We are currently interested in the heavy-light mesons described by the field H- . 
( It should also be possible to get static solutions, solitons, that deviate from the vacuum 
by a finite amount, describing for example baryons with one heavy quark). 

Expanding around the vacuum, we find that the equation for H decouples from the 
others to linear order: 



.dHHry) 1„ 

% = —a 

du 2 



J ds'Hi(r,s')G(s , ,r)5_(s , ,r , ) + J ds' Hi (r, s')h(s' ,r'). (33) 
The constraint on H becomes to linear order, 

J H!(r,r')6+(r',r")dr' = 0. (34) 
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Now, H(r',r) is the wavefunction of the heavy-light meson; we should expect this to 
be the product of a heavy quark wavefunction concentrated at the origin, and a light quark 
wave function. The ansatz H-(r\r) = 8(r f )ip(r)c! is therefore reasonable. The constant 
c\ will determine the internal quantum numbers ( spin,flavor) that label the degenerate 
levels of the meson. We get 

= ^aj ds'i;(s')G(0,s')S-(s',r') + J ds'^{s')h{s' ,r'). (35) 

Now recall that h(s',r') = —h(r',s'), 5-(s',r') = 8 + (r',s') so that 



dib(r') 1 „ 

—i — = — a 

du 2 



ds'S + (r',s')i>(s')G(s',0) + J ds'h(r',s')i>(s'). (36) 



Thus we arrive at the equation for stationary states, 

~[id r + k 2 r- 1 (id r )- 1 r- 1 ]il;(r) + J 8+(r, s)V(s)ip(s)ds = Eip. (37) 

Here, V(r) = — ^ is the Coulomb potential and E the binding energy of the meson. 
This is to be supplemented by the constraint 

8-(r,s)ip(s)ds = 0. (38) 

The above pair of integral equations will determine the meson wavefunction and 
masses in our model. They describe the propagation of a light quark in the spherically 
symmetric potential created by the heavy quark. Since there are no gluon self-interactions 
within our approximations, this potential has turned out to be the Coulomb potential. 
Therefore the eigenvalue problem we get is scale invariant. 

In fact we should expect the scale invariance to be broken by quantum effects; no 
scale invariant equation can have a discrete spectrum of bound states. In the next section 
we will perform a renormalization of the coupling constant which will make sure that the 
ground state energy is finite. 

If we expand 

/■OO 

^(r) = / u x (r)^(X)d\ (39) 



o 



the constraint is automatically satisfied. The remaining equation can also be simplified, 

1 f 00 

A^(A) + -a / K(X, X')?P(X')dX' = E^{X) for A > 0. (40) 

2 Jo 
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Here 

K(X,X') = J drV(r)u* x (r)u y (r) (41) 

is the integral kernel of the potential in the basis diagonalizing h. This is the analogue of 't 
Hooft's integral equation. Using the explicit form of the functions u we can simplify these 
equations further. (For example if k = 1 they are spherical Bessel functions; for arbitrary 
k, u(r) can be obtained as a power series.) Our result is then similar to a Wiener-Hopf 
integral equation. The two dimensional analogue of this equation is discussed in ref. [13]. 

5. Renormalization and Beta Function 

The eigenvalue equation of the free hamiltonian 

~{id r + k 2 r- 1 (id r )- 1 r- 1 )u(\r) = ^u(Xr) (42) 

can also written as the differential equation 

k 2 

[-d r {rd r u) + —u] = -iXdJru). (43) 
r 

The solutions are of the form e %Xr times a polynomial in jj^. For the special case k = 1 
which we will now study in detail, 

= " (44) 

These functions are orthonormal in the sense that 

poo 

/ u*(Xr)u(X'r)dr = S(X — A'). (45) 
Jo 

The solutions to the constraints are wavefunctions of positive radial momentum: 

POO 

^(r) = / 4>(X)u(Xr)dX (46) 

The integral equation is equivalent to minimizing the the energy, 

, r 7l _ / oo A|^(A)| 2 dA + / oo y(r)|^(r)| 2 dr 

where 1 1^| | 2 = J °° |^(A) \ 2 dX. 

Since it is difficult to solve this problem exactly, we will find the ground state energy 
by a variational principle. 
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However there is a problem: the ground state energy is divergent. The potential and 
kinetic energies are both of dimension 1 in energy units, so that with ^(A) = ^(/^A) 

£(VV) = (48) 
A* 

for any \x > 0. If there is one state with E < ( an example of which is given below), 
there are states of arbitrarily negative energy. We must introduce a cutoff and make the 
bare coupling coupling constant depend on it such that the ground state energy remains 
finite as the cutoff is removed. If a > is such a short distance cutoff, The regularized 
energy will be of the form 

E(a, a (o)) = a~ 1 E 1 (a (a)) (49) 
by ordinary dimensional analysis. Define the beta function by 

P(& (a)) = -o-^- (50) 
The condition that E be independent of a is then determines the beta function: 

/3(ao) = -^)- (51) 

Such a nonperturbative renormalization scheme was introduced by Thorn in a nonrela- 
tivistic context [14]. 

As in field theory, a naive short distance cutoff is not convenient. It is better to use 
instead an analogue of dimensional ( or analytic ) regularization. Define the regularized 
potential to be 

v " (r) = ~T (52) 

We could then compute the regularized energy, 

E(no,a ,e) = (i Ei(a,e) (53) 

from which the beta function can be determined. 

We are not able to calculate the regularized ground state energy exactly; instead, we 
will estimate it by a variational principle. The ansatz 

^(A) = \ e - bX for A > (54) 
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gives 

E< I- -^{-+4 log mu b-3}+ 0(e) (55) 
lb no e 

Minimizing the r.h.s. in the variational parameter b gives 

W=(1 _I + |] + |_ (56, 

so that 

7rb(a ) 

This gives, 

« 5 > " -ra- (58) 

This beta function has only one zero at a = 0, near which 

M = _^A + ... (59 ) 

so that our renormalized theory is asymptotically free. This can be compared with the 
well known one-loop result of perturbation theory [15], [16], [17], [18] 

0(a) _ 11 iV 2 - 1 a 
a 3 2N vr 1 ' 

Recalling that a = Na we see that ( for large N) we have a 16 where the one- loop result 
has a 11; i.e., agreement to about 30 %. 

We found that a more general ansatz such as 

ip(X) = A(l + cA)e~ bA (61) 

has a minimum at c = 0; i.e., does not improve the energy. The ansatz 

^(A) = A c e" 6A (62) 

has a minimum at c = 1.17 and lowers the energy by less than a percent. To get a 
substantially better estimate of the ground state energy (and hence the beta function) we 
must solve the problem numerically. It is also of much interest to study the excited states. 
These issues are currently under study and we hope to report on them soon. 
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6. Appendix: Constraints on Fermion bilinears 

In this appendix we will describe how the large N c limit of fermion bilinears leads 
to a classical theory whose phase space is the Grassmannian. One can think of this as a 
quantum mechanical analogue of 'bosonization'. 

Let us consider a set of operators satisfying Canonical anticommutation relations 
(CAR): 

[x tao ,X/36]+ = W (63) 

all other pairs of anticommutators being zero. Here a, (3 = 1 • • • N c while a, b are all others 
quantum numbers such as spin, flavor, momentum etc. It will not matter to us for now 
what exactly they represent. For simplicity we will assume that they have a finite range. A 
representation for these CAR can be constructed on the Fermionic Fock space T as usual. 
Now consider the subspace JF of 'color singlet' states: 

llx ]aa ,Xpa]\s>=0 for | S >G^o- (64) 

(The operators ^[x^ aa ,X/3 a ] generate the algebra U(N C ) of 'color'.) In particular this 
condition implies that exactly half the one-particle states are occupied. 

The color singlet bilinears 

^[xKxab] (65) 
map Tq to itself. In this subspace, they obey the linear constraint 

K = 0. (66) 

Also, they satisfy the algebra 

The meaning of this is that N c &% form a representation of the Lie algebra of the Unitary 
group on JF . This is an irreducible representation so that the only operators in JF that 
commute with all the are multiples of the identity. 

The form a complete set of observables of the quantum system whose Hilbert space 
is Fq. Within the space of color singlet states, any observable can be written in terms of 
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the However, they are not independent variables; they satisfy a quadratic constraint 
between states in JF : 

= {\ + + 2^«- (68) 

To derive these one just has to reorder the factors so that the color singlet condition on 
the states can be used. 

In the large N c limit this becomes just 

= (69) 

Now, if we define the normal ordered bilinear 

we will get in the large N c limit the constraint 

M a c M c h - i(M c °eg + e«M|) = 0. (70) 

To derive the constraints in the text we must replace x by (^Q^j an< ^ the matrix <E> by 

(M-\t #t \ 
V H P-\S)' 

In the large iV c limit we get as dynamical variable a hermitian traceless matrix whose 
square is a multiple of the identity. The phase space, which is the set of all such matrices, 
is the 'Grassmannian'. The commutation relations now are replaced by Poisson brackets 
on the classical variables. We can recover the finite iV c theory by 'quantizing' this classical 
theory whose phase space is the Grassmannian. The wavefunctions of the quantum theory 
can be chosen to be some sort of functions on the Grassmannian. The precise statement 
is that the wavefunctions are holomorphic sections of a line bundle on the Grassmannian. ( 
The Grassmannian is a complex manifold, so the concept of holomorphicity makes sense.) 
This description of the wavefunctions is analogous to the coherent state picture of the 
harmonic oscillator. Line bundles on the Grassmannian are labelled by an integer ( Chern 
class) which can be identified with N c . For each N c therefore we have one quantum theory; 
the Hilbert space of this quantum theory can be shown to be just JF . This way we can 
recover the finite N c theory as a quantization of the classical theory on the Grassmannian. 
In particular this construction shows that there are no other constraints on the in the 
large N c limit. 
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We will encounter some divergences in extending this to the infinite dimensional case. 
However once we deal with normal ordered operators, the divergences can be handled 
consistently. A more formal approach is to define the infinite dimensional Grassmannian 
with certain convergence conditions and to construct the space JF as a space of sections 
of line bundles on it. [19] [20]. We do not need this construction for the purposes of this 
paper. 
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